data <- read.csv("./RichardsonRK4.csv",head=FALSE,sep=",")
# now can access each as data$V1 or data$V2

pdf("RichardsonOrbits.pdf") # makes a pdf printing device for plots
plot(data$V1,data$V2, pch=".",ylab="Position",xlab="Time",cex=0.8)
dev.off() #closes the pdf printing device

print("Updated ./RichardsonOrbits.pdf")